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ABSTRACT 

Although white dwarfs are believed to be the end point of most stellar evolution, unlike main sequence stars, 
they have not yet been the subject of dedicated time-domain surveys for exoplanets. We discuss how their 
size and distinctive colour make them excellent targets for wide-field searches for exoplanets. In particular, we 
note that planets of Earth-size can give rise to multi-magnitude eclipses of massive white dwarfs. Such a large 
signal is almost unmistakable and would be detectable even with very low-precision photometry. For objects 
of smaller size, the high accuracy photometry currently being used to detect Super-Earth and smaller planets 
transiting Rq stars, is capable of revealing minor planets down to R ~ 100 km as they transit white dwarfs. 
Such observations can be used to test current evidence for asteroid-size objects being the cause for dust rings 
which have recently been observed for a number of white dwarfs. No other current exoplanet search method is 
capable of detecting such exo-asteroids. 

As an initial test of this search strategy, we combine synoptic data from the Catalina Sky Survey with multi- 
colour photometry and spectra from the Sloan Digital Sky Survey to search ~ 12,000 white dwarf lightcurves 
for eclipsing events. We find 20 new eclipsing white dwarf binary systems with low-mass companions. This 
doubles the number of known eclipsing white dwarfs and is expected to enable the determination of accurate 
white dwarf radii. Three of the discoveries have radii consistent with substellar systems and show no evidence 
of flux from the eclipsing object in their SDSS optical spectra, or near-IR data. 

We discuss how future deep wide-field surveys, such as LSST, should reveal thousands of eclipsing white 
dwarf systems and thus put very strong limits on the asymptotic giant branch (AGB) survival of terrestrial 
planets around the progenitors of the present Galactic white dwarf population. 

Subject headings: general — stars: general — stars: binaries:eclipsing — stars: white dwarfs — planetary 
systems — minor planets, asteroids 



1. INTRODUCTION 

Since the discovery of planets with masses similar to Earth 
in orbit of Pulsar PSR 1257+12 (Wolszczan & Frail 1992; 
Wolszczan 1994), groups of astronomers have been attempt- 
ing to discover similar objects in the more habitable environs 
of main-sequence stars. Efforts towards this end were ini- 
tially largely concentrated on finding the very small radial 
velocity variations as a planet orbits its parent star (Mayor 
& Queloz 1995; Butler & Marcey 1996). The amount of 
time dedicated to such surveys remains reflected in the fact 
that the largest number of planets have been discovered in 
this way0 However, although the information from such sur- 
veys has undoubtedly increased interest in planetary searches 
and led to constraints on the distribution of massive planets, 
until now such surveys have required significant time with 
large telescopes and generally only leads to mass, period and 
semi-major axis constraints on individual planets (Patel et al. 
2007). 

For the past decade many groups have begun to search for 
planets around stars by looking for very small dips in bright- 
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ness as a planet transits its star. The presence of transits along 
with radial velocity enables the determination of density, and 
further observations during transits can constrain the atmo- 
spheric composition of the transiting object (Swain, Vasisht, 
& Tinetti 2008; Sing & Lopez-Morales 2009). Surveys have 
also been been successful in finding the microlensing sig- 
nature of exoplanets (e.g., Beaulieu et al. 2006). However, 
because of the accuracy required to detect Earth-size planets 
with transits, microlensing or radial velocities, such discover- 
ies remain elusive. The searches for low-mass planets in tran- 
sit surveys are now turning to the extreme accuracy possible 
with space-based observatories, such as Kepler (Batalha et al. 
2002) and CoRoT (Rouan et al. 1998). Ground-based transit 
searches are beginning to turn to M-dwarfs as exoplanet tar- 
gets, since a transit of an M-dwarf with radius 0.25 Rq gives 
rise to a signal 16 times that of a IRq star (Shankland et al. 
2006). However, even in this case the detection of Earth-sized 
planets is challenging for ground-based photometric searches 
(Montgomery & Laughlin 2009) 

1.1. White Dwarfs as Planet Hosts 

White dwarfs have radii 10 to 20 times smaller than M- 
dwarfs. They are the end point in the evolution of most stars 
(M < 8M Q ) and can provide insight into the history of star 
formation in the Galaxy as well as the fate of our solar system. 

It has long been theorized that close companions to the pre- 
cursors of white dwarfs do not survive the passage through the 
asymptotic giant branch (AGB) as objects within ~ 5 AU un- 
dergo a common envelope phase where the material is trans- 
ferred from one object to the other. Such systems are believed 
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to lead to cataclysmic variable (CV) systems where the com- 
panion star accretes matter onto a disk surrounding the white 
dwarf. In cases where the companion lies well beyond the 
limit of the expansion of the white dwarf precursor the or- 
bital distance expands as the system radiates energy (Livio & 
Soker 1984) and the expanded star loses mass (e.g., Schroder 
& Connon Smith 2008). Villaver & Livio (2007) theorized 
that gas planets smaller than Mj will not survive when they 
are within orbital distances of 3-5 AU of a main sequence 
star evolving through the AGB and planetary nebular phases. 
While, on the other hand, gas planets more massive than 2 Mj 
can survive down to ~ 3 AU. In some cases, planets may ac- 
crete mass and become a low-mass companion. However, as 
Villaver & Livio (2007) note, the maximum mass of a planet 
that can be evaporated inside the envelope of an AGB star de- 
pends on the efficiency with which the envelope is ejected, 
and this is not known. Clearly therefore, the limits given 
therein are very uncertain. In contrast to Villaver & Livio's 
(2007) result, Duncan & Lissaur (1998) found that the terres- 
trial planets of our planetary system (other than Mercury) may 
survive the AGB phase in orbits similar to their current ones 
(moved adiabatically outward by an amount depending on the 
Solar mass loss; see also Schroder & Connon Smith 2008). 
Indeed, Rasio (1996) also demonstrated the importance of the 
tidal dissipation in theories where objects such as Earth do not 
survive the AGB phase. 

In contrast, Debes & Sigurdsson (2002) noted that although 
planets within 5AU are expected to be swept up in the AGB 
phase, the significant mass loss that occurs will cause plan- 
ets to move in their orbits. In systems with multiple planets, 
they found that the interaction between planets orbiting a mass 
losing stars will cause some of the orbits to become unstable. 
This theory is similar to a proposed explanation for the large 
number of hot-jupiters discovered in small orbits around main 
sequence stars. Studies of the interaction in two Jupiter sys- 
tems by Ford, Havlickova & Rasio (2001) found that for pairs 
of planets beyond 5 AU, a large fraction would settle into a 
configuration where planets were significantly closer than the 
initial system. Jura (2008) considered the destruction of mi- 
nor planets during evolution of stars with a ranges of masses 
and predicted orbital radii at which rocky objects of various 
sizes are expected to survive to the white dwarf evolutionary 
phase. Clearly, there is still significant uncertainty in whether 
some rocky planets can survive the evolution of their host star. 

Another path by which a planet might survive or side- 
step AGB evolution is in the form of an sdB star. Bear & 
Soker (2010) note that blue subdwarf (sdB)/extreme horizon- 
tal branch (EHB) stars are prime targets for planets. Since 
EHB/sdB/sdO stars evolve to AGB-manque stars (Greggio 
& Renzini 1990) or post-early AGB (Brocato et al. 1990) 
stars, these stars do not evolve to the AGB stage before evolv- 
ing to the WD cooling track (see, e.g., Dorman, Rood, & 
O'Connell 1993 and Catelan 2007 for reviews on the sub- 
ject). All such stars with planet systems are thus expected to 
inevitably evolve into WD plus planet systems. Evidence for 
this evolutionary path was obtained by Silvotti et al. (2007) 
who discovered a 3.2Mj planet orbiting sdB star and note this 
as evidence that a planets can survive to orbit a post red gi- 
ant star at distances less than 2AU. Additional similar sys- 
tems have been discovered by Geier et al. (2009) and Lee et 
al. (2009). The resulting WDs of these systems should have 
relatively low-mass planets orbiting them. In like vein, if the 
red giant branch (RGB) progenitor of the sdB stars were to 
lose additional mass during the RGB ascent, they would be- 



come He-core WDs (Catelan 2007, and references therein). 
The resulting WDs would clearly have relatively low-mass 
planets orbiting them. 

Apart from surviving the AGB, close companions to WDs 
should also survive the RGB ascent of the progenitor star. In- 
deed, the fact that there exist planetary systems around sdB 
stars (e.g., Silvotti et al. 2007; Geier et al. 2009; Lee et al. 
2009) conclusively demonstrates that many planets are able to 
survive the first ascent of the giant branch by their host stars. 
During such an ascent of the giant branch, the envelope ex- 
pands to tens of solar radii and could consume a very near 
companion. The recent theoretical work of Villaver & Livio 
(2009) predicts that planets can survive the RGB evolution 
while residing in close proximity to the stellar host in the cor- 
rect mass range (see their figure 1). They show that low-mass 
companions are favored to survive near the stars. Villaver & 
Livio (2009) predict that a Jupiter-mass planet can avoid cap- 
ture by their parent star within a 0. 18 AU orbit - a range that 
is suitable for WD progenitors. 

1.2. Evidence of RGB and AGB Survival 

Currently surveys for exoplanets have discovered more than 
350 planets. As the number of stars expected to have a gas 
giant planet increases linearly with mass in the 0.4- 3M Q 
range (Kennedy & Kenyon 2008), it is not surprising that 
Jupiters are commonly found orbiting the progenitors to WDs. 
Indeed, planets have been discovered around stars evolving 
along the RGB. For example, de Medeiros et al. (2009) dis- 
covered a 1.9 - 2 Mj object orbiting the 2M Kill giant star 
HD 110014. Additional discoveries of planetary compan- 
ions to evolved stars have been made by Peek et al. (2009), 
Valenti et al. (2009), and Johnson et al. (2007, 2008). Bower 
et al. (2009) found that 26% of evolved A-type stars should 
have Jovian companions within 3 AU. 

Although stellar evolutionary theory does not favor the sur- 
vival of planets in close proximity to a star evolving through 
the AGB phase, there is mounting evidence of low-mass 
stellar companions to white dwarfs. For instance, Fahiri, 
Burleigh, & Hoard (2008) recently discovered an L8 dwarf 
star orbiting a white dwarf, and Unda-Sanzana et al. (2008) 
found a similar system. In these cases, the companions are 
believed to be in interacting CV systems. The detection of 
substellar companions to CVs is complicated by the fact that 
such systems are usually detected via variability caused by 
accretion processes. This effectively rules out the possibil- 
ity of determining the original mass of the companion star. 
Nevertheless, it has been suggested that the secondaries in 
CV systems will end up as planetary mass objects (Patterson 
1998) or low-mass L dwarf-like stars (Howell, Rappaport, & 
Politano 1997). In addition, Kilic et al. (2007) discovered a 
companion to a white dwarf with mass < 0.1 M Q , ruling out 
a main-sequence companion. Burleigh et al. (2006) found an 
L8 brown dwarf companion to white dwarf WD0137-349B, 
and Steele et al.(2009) discovered a similar L dwarf orbit- 
ing white dwarf PHL-5038A. Based on arrival time arguments 
and Spitzer data, Mullally et al. (2009) have recently placed 
an upper limit on the mass of the substellar companion to an- 
other white dwarf, GD-66, as 5-7Mj. 

Evidence for the material comprising terrestrial planets sur- 
viving the AGB evolutionary phase comes from the discov- 
ery of white dwarfs with dusty disks. These disks have been 
attributed to recent disruption of minor planets and comets 
(Mullally et al. 2007, and references therein), as metals from 
the accreted asteroids should only be seen in white dwarf 
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spectra for a few millions of years before settling to the core. 
Zuckerman et al. (2007) discovered heavy elements in the 
spectrum of white dwarf GD 362 and proposed it to be the 
remnant of an asteroid with a composition similar to Earth's 
crust that "may have once been part of a larger parent body 
not unlike one of the terrestrial planets of our solar system." 
Brinkworth et al. (2009) also discovered a dusty component 
to the gaseous disk orbiting white dwarf SDSS J1228+1040. 
Koester (2009) calculated the diffusion timescales for various 
types of elements accreted on white dwarfs and found that, 
for GD 362, the heavy elements would gravitationally settle in 
less than 10 6 years. Since the lifetime of a debris disk from a 
disrputed asteroid is < 10 6 yr, the observed disruption must be 
a recent event that took place long after the AGB phase (Jura 
2008; Kilic et al. 2008). For this to occur the sources of rocky 
material must have survived the AGB as well as the Planetary 
Nebula evolutionary phases before their final encounter with 
the white dwarf. 

Recently, Klein et al. (2010) discovered evidence for accre- 
tion of rocky material from an extrasolar minor planet in the 
spectroscopic analysis of the atmosphere of He-atmosphere 
white dwarf GD 40. Melis et al. (2010) discovered evidence 
for dusty and gaseous disks around three white dwarfs which 
they attributed to the remnants of planetary systems that sur- 
vived AGB evolution. Similarly, Dufour et al. (2010) dis- 
covered a He-atmosphere white dwarf with pollutant mate- 
rials having relative abundances resembling the composition 
of Earth as well as a circumstellar disk. The authors suggested 
that the amount of material was due to accretion from an ex- 
trasolar body that was as large as Ceres. Hogan et al. (2009) 
note that the existence of metal-rich dust disks suggests that 
even terrestrial planets can survive the final stages of stellar 
evolution. Furthermore, Jura et al. (2009b) used XMM and 
Spitzer to monitor G29-38 and GD 362. G29-38 was discov- 
ered by Zuckerman & Becklin (1987a) and originally believed 
to contain a brown dwarf companion. Jura et al. (2009b) 
found it to be accreting heavy elements from a low-mass body 
and that GD 362 is accreting material from either a large num- 
ber of Ceres-like asteroids or a single body that was probably 
as massive as Mars. Six more white dwarfs with dusty silicate 
rings were found by Jura, Farihi, & Zuckerman (2009), thus 
providing additional proof that planetary material can survive 
AGB evolution and end up in orbit of a white dwarf. Lastly, 
as noted above, when a white dwarf evolves from an sdB 
star, it is unnecessary for planetary companions to survive the 
AGB phase of the star, since the AGB phase is completely 
bypassed. 

1.3. Searches for Planets around White Dwarfs 

Many searches have been carried out to discover low-mass 
companions to white dwarfs. These have largely involved 
searches sensitive to the near-IR light, as white dwarfs, in con- 
trast to Jupiter-like planets, emit little light at infra-red wave- 
lengths. Probst (1983a, 1983b) and Zuckerman & Becklin 
(1987a, 1987b) searched for cool companions to 100 WDs 
using near-IR J,H,K data. Zuckerman & Becklin (1988, 
1992) discovered the first candidate brown dwarf orbiting 
white dwarf GD 165. Wachter et al. (2003) and Hoard et 
al. (2007) searched for late-M- to early-L-dwarf companions 
using 2MASS sky survey data, but found very few objects 
had potential brown dwarf companions. Farihi et al. (2005a) 
carried out a proper-motion plus near-IR survey and discov- 
ered a second brown dwarf, GD 1400 (Farihi et al. 2005b), 



with an implied orbital separation of < 1 AU (Burleigh et 
al. 2006). Steele et al. (2009) used H- and K-grism spectra of 
PHL 5038 and discovered the third brown dwarf companion 
known to orbit a WD. As part of the DODO survey, Hogan, 
Burleigh, & Clarke (2009) used J-band observations to search 
for very low-mass brown dwarfs and giant extrasolar plan- 
ets as common proper motion companions to 23 nearby white 
dwarfs. No planet counterparts were found within orbital sep- 
arations in the range of 60 - 200 AU. Debes et al. (2005, 
2006) searched for sub-stellar companions to a total of 20 
WDs using HST NICMOS coronagraphic images, CFHT J- 
and H-band data. An additional brown dwarf companion was 
found in the near-IR spectrum of WD 0137-349 by Burleigh 
et al. (2006). This object has a period of less than 2 hours, 
suggesting that it had either survived common envelope evo- 
lution or migrated to a short orbital period. 

Radiszcz (2008) searched for low-mass brown dwarf and 
massive giant planets in VLT/NACO and Magellan/PANIC 
data for 116 nearby white dwarfs, and concluded that substel- 
lar objects more massive than 10Mj are rare, with an upper 
limit of < 5% of his sample stars harboring any such com- 
panions. Debes et al. (2007) obtained Spitzer IRAC data 
for 5 DA WDs in order to search for mid-IR excesses. Simi- 
larly, Mullally et al. (2007) used Spitzer to look for planetary 
companions and disks around 124 white dwarfs and found 
one object consistent with a brown dwarf companion. Farihi, 
Becklin, & Zuckerman (2008a) also examined Spitzer IRAC 
data for 48 WDs and discovered no new companions, con- 
cluding that current evidence suggests that white dwarf bi- 
nary systems with L-type companions are rare - although, as 
noted earlier, some such systems have indeed been found. For 
six of the WDs this survey was sensitive to massive plane- 
tary companions with masses > 5Mj, and for 20 WDs it was 
sensitive to planetary companions with masses > 10Mj. The 
best evidence for a planetary mass companion was discov- 
ered by Mullally et al. (2008), who used timing experiments 
to search for period variations among 15 known variable DA 
white dwarfs. They found evidence for a 2Mj planet in a 4.5- 
yr orbit of GD 66, corresponding to a 2.3 AU orbit. 

Gould & Kilic (2008) noted that current planet searches do 
not probe planets around massive stars and that such searches 
could be extended to higher masses by searching for planets 
around the bright young white dwarfs that are the remnants 
of massive stars (3M Q < M < 7M0). Their proposed survey 
method was to once again use a Spitzer IR search to monitor 
49 young, hot, white dwarfs, a very small fraction of the many 
thousands of known white dwarfs in the online version of the 
McCook & Sion (1999; hereafter MS99) catalog]] Kilic et 
al. (2010a) noted that the the lack of planetary discoveries 
may have been due to the fact that WDs in randomly targeted 
searches are mainly the remnants of low-metalicity progeni- 
tors. 

Other planetary search methods are used to find very close 
planetary companions to much brighter stars. These include 
radial velocity, transit, and eclipse searches. The broad ab- 
sorption lines of white dwarfs are not ideal for obtaining ac- 
curate radial velocities, and since white dwarfs are gener- 
ally faint, such measurements would be very time consum- 
ing. Photometric searches for transits and eclipses can survey 
much fainter target stars more efficiently. 

1.4. Detecting Planetary Eclipses and Transits 
7 http ://www. astronomy. villanova. edu/WDCatalog/index .html 
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The exact sizes of white dwarfs are not yet very well con- 
strained by empirical results (Pyrzas et al. 2009), but those 
with measured sizes vary from slightly smaller than Earth for 
the most massive, to a couple of times Earth's size for the 
least massive. A summary of the 14 known white dwarf main 
sequence eclipsing binaries is given by Pyrzas et al. (2009, 
their Table 9). In short, these binaries have orbital periods be- 
tween one half and one eighth of a day and secondaries that 
range from spectral types K2 to M6.5. The effective temper- 
atures of these white dwarfs range from 7,000 to 57,000 K, 
their masses range from 0.4 to 0.8M Q , and their radii are in 
the range 0.009- 0.022 R Q . 

The diminutive nature of white dwarfs means that unlike 
main sequence stars, a planet as small as Earth can com- 
pletely eclipse a white dwarf. Unlike transit surveys which 
require milli-magnitude precision to discover transiting plan- 
ets, the signal of a planetary eclipse amounts to many magni- 
tudes, and increases with decreasing luminosity of the eclips- 
ing object. A complete eclipse caused by an Earth-like planet 
would cause the white dwarf to drop beyond the limits of de- 
tectability by even the largest current telescopes. The detec- 
tion of such a signal could thus be possible even with very 
poor-quality photometry. Furthermore, objects the size of 
the Moon (0.27 R ffi ) would cause short, but easily detectable, 
planetary transits with depth equivalent to an M-dwarf eclips- 
ing a solar-type star. Although there has been increasing ef- 
fort concentrated on the discovery of transiting exoplanets, 
only one survey has been carried out to detect planets eclips- 
ing white dwarfs (Wilson 2003). In that survey, only 75 white 
dwarfs were followed. 

For the purposes of comparison, we note that the recently 
discovered exoplanet COROT-Exo-7b has size ~ 1.7R© (or 
0.015^0) and orbits a G9V star with a radius ~ 0.85/?© 
(Leger et al. 2009). The white dwarf GD 362 has mass 
0.73 M© and is expected to have a radius of ~ 0.01 1R©. 
Therefore, an object with a radius of 140 km orbiting a white 
dwarf like GD 362 would cause a transit signal of equivalent 
depth to COROT-Exo-7b transits. The source of GD 362's 
dust ring is proposed to be the remains of a recently disrupted 
Mars-size object, or many Ceres-size asteroids (Jura et al. 
2009b). Indeed, Koester (2009) estimated that the amount 
of material accreted on GD 362 over the lifetime of a de- 
bris disk (~ 10 5 years) would be similar to that of the dwarf 
planet Ceres (radius ~ 450 km). Therefore, it is possible to 
detect minor planets orbiting nearby white dwarfs with mis- 
sions such as Kepler. The large signal of transits or eclipses 
of terrestrial planets around WDs means that a much larger 
volume can be surveyed than with main-sequence stars. 

Assuming Jupiter-like companions to white dwarfs occur 
with the same distribution as in solar-like stars, the line-of- 
sight requirements for eclipses means that Jupiter-size planet 
eclipses of white dwarfs would be 10 times less common than 
Jupiter-like transits of solar-sized stars. Similarly, Earth-size 
planet eclipses would be 100 times less common than tran- 
sits of solar-type stars. In Figure Q] we present the proba- 
bility of observing an eclipse along any given line-of-sight 
for objects from Earth-size to Jupiter-size eclipsing an Earth- 
sized white dwarf (similar to Sirius-B) using a simple geo- 
metric model. Very close orbits with high eclipse probabili- 
ties can be ruled out as being unstable within the Roche limit 
(although dependent on composition). For planets similar to 
solar system ones, the Roche limit varies from ~ 0.004AU 
for a rigid planet of Earth's density to 0.011 AU for a fluid 
planet the density of Jupiter. Of the 56 transiting planets 



with currently known orbital radii, 68% have orbital radii less 
than 0.05 AU. For Jupiter-size planets orbiting a white dwarf 
within this range, there would be > 1% chance that it would 
be observed to eclipse. Nevertheless, Figure 1 demonstrates 
the need to monitor thousands of white dwarfs to observe 
eclipses of Earth-sized planets. 

In Figure [2] we show the relationship between orbital pe- 
riod and eclipse length for objects ranging from late-M dwarfs 
with radii of 0.25/?© and mass 0.2M© to Earth-like objects. 
Clearly, objects with eclipses lasting less than 20 minutes are 
the best substellar candidates. 

In this paper we first discuss the collection of data from the 
SDSS (York et al. 2000), 2dF (Croom et al. 2004) and MS99 
catalogs to search for candidate white dwarf stars. Next, we 
discuss matching these objects to Catalina Sky Survey syn- 
optic data where light curves of these objects are built. We 
then discuss the determination of the periods of the eclipsing 
objects and place limits on the companions based on near-IR 
photometry, the phased light curves and spectra. Finally, we 
summarize the results. 

2. WHITE DWARF SELECTION 

As demonstrated in the previous section, to discover white 
dwarfs being eclipsed by their companions, it is necessary 
to have both a large sample of white dwarfs and significant 
amount of data covering each object. The Catalina Sky Sur- 
vey (hereafter CSS) has used the 0.7m f/1.9 Catalina Schmidt 
Telescope to discover Near-Earth Objects (NEOs) since mid- 
2005. The eight square degree field of view of the Catalina 4k 
x 4k camera makes it ideal for the discovery of optical tran- 
sients of all kinds (Drake et al. 2009). On a clear night, the 
Catalina Schmidt Telescope typically covers ~ 1200 square 
degrees of sky in a sequence of four 30-second exposures 
from set fields covering > 30,000 square degrees in the dec- 
lination range -40° < 5 < 70°. Observations are unfiltered 
and generally limited to Galactic latitudes \b\ > 10° to pre- 
vent confusion caused by crowding. Weather permitting, ob- 
servations have been made 21 nights per month, and typi- 
cally reach V magnitudes of 19 to 20. Four images of a field 
are taken in a sequence spaced evenly over approximately 30 
minutes. The CSS dataset consists of fields covered from a 
few times to more than 250 times. Since white dwarf eclipse 
timescales, Te, are 10-30 minutes, CSS data are excellent to 
detect eclipsed white dwarfs. However, the lack of any colour 
or spectroscopic information from CSS makes it necessary 
to select white dwarfs using other large surveys, such as the 
SDSS. 

Eisenstein et al. (2006) recently produced a catalog of 
~ 9000 spectroscopically identified white dwarfs from the 
SDSS DR4 data release. For the first part of our analysis 
we matched the Eisenstein et al. white dwarf catalog to ob- 
jects covered by CSS. The SDSS catalog reaches white dwarfs 
as faint as r = 21.4. Since we need to discover significant 
eclipses in the CSS photometry, which reaches CSS r ~ 20, 
we selected only matching SDSS objects with magnitudes 
u < 18.5 or r < 18.6, reducing the number to 2684 white 
dwarfs, or ~ 29% of the Eisenstein et al. catalog. Among 
these objects we further culled the list by selecting only ob- 
jects covered by CSS data on more than five nights (> 20 im- 
ages). This list consisted of 1258 white dwarfs. For each 
white dwarf, we produced a light curve and a series of im- 
age cutouts. We then examined each light curve and checked 
images associated with outliers. Many of the white dwarfs ap- 
peared to be variable at some level. Some of these are likely 
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to be ZZ Ceti, DBV, or DOV variables, while a few others 
may be AM CVn variables. It was clear that only about 10% 
of the light curves exhibit potential eclipses, and almost all of 
these were due to image artifacts, blends, image edge effects 
and bad pixels. 

We also included the 227 white dwarfs spectroscopically 
identified by the 2dF project (Croom et al. 2004). Only 45 
of these fall in the well-covered CSS fields. Additionally, 
the 112 cool white dwarfs discovered by Kilic et al. (2006) 
were examined. Of these, 89 were covered by CSS. Further- 
more, we selected white dwarfs from the MS99 database and 
removed those already included from the SDSS catalog. As 
most of the remaining MS99 objects have very poorly con- 
strained positions, we matched these objects with the USNO- 
B catalog and selected the nearest blue star to the MS99 cata- 
log location. 1843 white dwarfs were selected from the MS99 
database after removing the SDSS matches. We selected stars 
with magnitudes 10 < m < 19 and B-R < 0.4 in either set 
of image epochs within a 22" by 14" box. 567 of these were 
well covered by CSS data. The total number of confirmed 
WD covered by CSS data that we checked was thus 1959. 

2.1. Colour Selection 

The second part of our white dwarf selection process in- 
volved searching for stars with colours consistent with white 
dwarfs using the SDSS data release six (DR6). The DR6 pho- 
tometry set (Adelman-McCarthy et al. 2008) covers more 
than twice the area of the SDSS-DR4 spectroscopic release 
used by Eisenstein et al. (2006). Our primary colour selection 
method among the SDSS-DR6 data is based on the colour- 
colour locations of white dwarfs confirmed by Eisenstein et 
al. (2006). Firstly, we selected objects with u-g < 0.6 and 
g — r < 0.2 and next, following the known white dwarf loca- 
tions, we cut out a region around the white dwarf sequence as 
shown in Figure [3] As our aim was to find white dwarfs with 
low-luminosity companions, rather than those with compan- 
ions that are bright enough to affect their observed colours, 
we were not concerned about removing objects with colours 
skewed by the presence of a luminous companion. After re- 
moving the objects included among the known DR4 white 
dwarfs, 18,211 candidate white dwarfs remained in this re- 
gion. We expect this colour selection will include many QSOs 
that overlap with white dwarfs in colour-colour space. How- 
ever, since QSOs are unlikely to exhibit eclipses, we expect 
little contamination among our candidates. Similarly, we ex- 
pect a number of sdB stars in this colour region. In Figure 
[4] we show the coordinates of all the white dwarf candidates. 
Of these, 12,994 had more than five nights of CSS observa- 
tions. Based on our initial test with spectroscopically identi- 
fied white dwarfs, we selected eclipsing candidates exhibiting 
at least one point 5a below the average magnitude, or with 
three points lying 3a below the average. This selection re- 
duced the number of candidate eclipses by a factor of ten. 
As with the spectroscopically confirmed white dwarfs, each 
of the images associated with outlier points was inspected to 
remove spurious candidates. 

3. RESULTS AND ANALYSIS 

In Table 1, we present 26 eclipsing white dwarfs found in 
this work, including 16 new discoveries. For each object, we 
also present the photometric measurements from the SDSS 
DR6 release and the near-IR magnitudes for those white 
dwarfs with matches in either 2MASS (Nikolaev et al. 2000) 
or UKIDSS (Lawrence et al. 2007). Proper motions from the 



USNO-B catalog (Monet et al. 2003) are also included. All 
the eclipsing systems, except for two, had spectroscopic con- 
firmation as WDs in public SDSS data or in earlier work. Five 
of the systems discovered are already known from work by 
Drake et al. (2009), Pyrzas et al. (2009), and Nebot Gomez - 
Moran et al. (2009). An additional four of the systems are 
known CV systems which can be clearly distinguished from 
non-accreting systems by the presence of emission lines in 
their spectra. 

In Figures [5] and [6] we present the location of the eclipsing 
white dwarfs in colour-colour space in relation to the spec- 
troscopically confirmed SDSS white dwarfs from Eisenstein 
et al. (2006). The presence of an M-dwarf secondary clearly 
influences the SDSS r—i and i — z colours. The few objects 
that lie on the main locus of crosses are the best candidates 
for systems with substellar companions. The groups of black 
points without overlapping crosses in Figure [6] are QSOs that 
lie within our initial white dwarf u- g vs g — r selection. This 
gives us an estimate of the number of QSOs falling within our 
colour selection. 

In Figures [7] and [8] we present the locations of the eclips- 
ing white dwarfs in relation to known white dwarf/M dwarf 
(WD-MD) pairs and CV systems. In Figure [7] we see that 
many WD-MD systems and CVs fall within our colour selec- 
tion. However, from Figure [8] it is clear that the r—i distri- 
bution is much closer to the pure white dwarf sequence than 
the WD-MD systems found spectroscopically by Silvestri et 
al. (2006). Not surprisingly, this suggests that Silvestri et al.'s 
spectroscopic selection may have missed many of the WD- 
MD systems where the M dwarf contributes little flux to the 
spectrum. Also it is clear that CV systems with disks and 
Balmer lines in emission are quite well separated from WD- 
MD pairs in g — r vs r—i colours. In Figure [8] we see that 
very few of the Silvestri et al. WD-MD systems have colours 
similar to those eclipsing systems discovered. This is further 
evidence that Silvestri et al.'s spectroscopic selection missed 
systems with faint M-dwarf companions. Since the proba- 
bility that a large M dwarf eclipses a white dwarf is much 
higher than that for a small, faint M dwarf, this suggests that 
most of the companions to white dwarfs are late-type objects, 
and most of these have been missed by Silvestri et al. be- 
cause their presence does not significantly affect the observed 
colours or spectra. 

During our search of the CSS lightcurves, we discovered 
a large number of slowly varying sources among the colour- 
selected white dwarf candidates. We examined SDSS spec- 
troscopic data for a number of these objects, and as expected, 
they were all QSOs. Objects exhibiting significant variability 
which was not caused by eclipses were noted separately from 
eclipsing objects. In Figure [9] we show the distribution of 
these objects in the (g-r) vs. (u-g) colour plane. Although 
most are probably QSOs and none show signs of outbursts, 
some of these are likely to be CVs which are varying signif- 
icantly due to accretion. Additional variable sources may be 
pulsating ZZ Ceti, DBV or DOV white dwarf variables. The 
QSOs remaining in our colour selection are a significant con- 
taminant to the determination of the total number of white 
dwarfs surveyed and the fraction of white dwarfs with com- 
panions. Of the ^5000 lightcurves inspected visually, 600 
exhibited light curve variations consistent with QSOs. Thus, 
we expect at least ~ 12% of the photometrically selected can- 
didates are QSOs. This number could be significantly larger 
if many of the QSOs exhibit very little variability. However, 
based on the overlapping fraction of variable objects in our 
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colour selected region, we expect less than 25% of the objects 
are QSOs. Another type of contaminant to our colour-based 
white dwarf selection are sdB stars. Based on the colours 
of these subdwarfs and the number identified by Kleinman et 
al. (2004) in SDSS data, we expect 10% of candidates are 
such objects. Correcting the total number of white dwarfs for 
the presence of QSOs and subdwarfs in our colour selection, 
we expect our total sample contains roughly 12,000 white 
dwarfs. 

3.1. Period Determination 

Period-finding techniques such as Lomb-Scargle peri- 
odogram analysis (Scargle 1982) or the analysis of variance 
(Schwarzenberg-Czerny 1989) do not work well when there 
are only a few points with significant signal, as happens in 
the case of the short eclipses of white dwarfs. To determine 
the periods for each of the eclipsing white dwarf systems, we 
used the simple method employed by Drake & Cook (2004). 
Here, as with phase-dispersion minimization (PDM, Stelling- 
werf 1978), the points are phase-folded and binned for a range 
of frequencies. However, unlike PDM, where the determina- 
tion can be biased by artificial trends within the data, only the 
points with significant signal (the eclipsed points) are folded 
until they lie within a single phase bin. Period aliases were 
removed by requiring that the eclipsed points form a continu- 
ous sequence in phase. That is, we enforce that the group of 
eclipsed points are not interrupted in phase by points near the 
baseline magnitude. 

Since the eclipse events last for only a fraction of the orbital 
period, it is necessary to have a large number of measure- 
ments to determine the periods of these objects. It is critical 
that the eclipses are measured on multiple nights spread over 
many days. For some systems there were insufficient data to 
uniquely determine the period of a system. 

As small companions will have eclipses shorter than larger 
objects for a given orbital period, there is a strong bias toward 
the larger objects being detected and having their periods de- 
termined. In Figure [TT] we present the phased light curves of 
two of the eclipsing white dwarfs, and in Table 2 we present 
the details of the eclipse periods determined from CSS pho- 
tometry. 

In FigureQ~2] we present the CSS and SDSS light curves of 
another eclipsing white dwarf. This white dwarf resides in the 
stripe 82 region that was observed more than 70 times by the 
SDSS-II supernova survey (Holtzman et al. 2008). This figure 
clearly demonstrates a modulation in brightness with phase 
that increases with wavelength. In the Gunn u-band the am- 
plitude is ~ 0.15, rising to ~ 0.35 in the z-band. Such modula- 
tions occur in binary systems when one of the stars is distorted 
into an ellipsoidal shape by the gravitational force of the other 
(Drake 2003). However, in this case the observed modula- 
tions occur once per orbit, rather than twice as expected for 
ellipsoidal modulations caused by a distorted secondary fill- 
ing its Roche lobe. Indeed, since the eclipsing object is more 
than two magnitudes fainter than the white dwarf, ellipsoidal 
modulations due to the secondary should be unnoticeable. In 
this case the most likely cause is the reflection effect, where 
the light from the white dwarf is reprocessed within the at- 
mosphere of the companion star and re-emitted. This effect 
explains the single modulation and the increasing effect with 
wavelength (Chen et al. 1995). 

The masses of most white dwarfs discovered by the SDSS 
range from ~ 0.3M Q to 1.1 Mq (Ciechanowska et al. 2007). 
As the SDSS spectroscopically confirmed white dwarfs fol- 



low the whole range of our colour selection, it is reasonable 
to assume that the eclipsing systems we discovered will also 
lie within this range. In Figure [TUl we present the eclipse 
times and periods for the systems discovered. The dashed 
line presents the eclipse times for a large O.3M white dwarf 
while the solid line shows the same quantity for a 1 . 1 M Q one. 
The figure suggests that most of the secondaries have sizes 
near 0.1 Rq- However, the times presented assume the sec- 
ondaries cause central eclipses. Near to the limb of eclipsing 
white dwarfs the eclipse time can be significantly shorter for 
companions that are much larger than the white dwarf. Nev- 
ertheless, a few objects exhibit very short eclipses consistent 
with very small companions. In FigureQI)] we also show dot- 
ted lines delimiting the CV period gap region. It has been well 
established that very few CV systems occur in the two to three 
hour period range (Davis et al. 2008). Few of the systems dis- 
covered appear in the CV period gap and thus, although most 
do not show the emission and disk signatures of CVs, there 
may be a related reason for this. Indeed, the companions may 
simply be objects too small to fill their Roche lobes until they 
are in very close proximity to the white dwarf. 

3.2. White Dwarf Distances 

In order to understand the nature of the eclipsing white 
dwarf systems, we determined the temperature and surface 
gravity of each of these white dwarf systems with an SDSS 
spectrum using the autofit program (Kleinman et al. 2004). 
Since the eclipsing objects we found are DA-type white 
dwarfs, next we apply the method of Holberg, Bergeron, & 
Gianninas (2008) to determine the spectral distances of each 
of the white dwarfs. Briefly, this involves interpolating the 
DA white dwarf models to determine the mass and radius of 
each object for a given T e s and log(g). We then correct the 
measured SDSS Gunn u- and g-band measurements for red- 
dening using the Schlegel, Finkbeiner, & Davis (1998) ex- 
tinction map, and use Holberg et al.'s (2008) SDSS absolute 
u- and g-band model magnitudes to determine the average dis- 
tance modulus of each object. We only use u and g photome- 
try, which span the wavelength range from 320 to 550nm, as 
the r, i and z photometry measurements could be biased by 
the presence of flux from a bright secondary. In Table 3, we 
present the white dwarf parameters determined from the spec- 
tra as well as the secondary types and radii determined from 
the CSS light curves and near-IR data, as described below. 

3.3. Companions 

We examined the spectra associated with each eclipsing 
white dwarf system to detect the presence of a secondary com- 
panion. In a number of cases the spectra show a very distinct 
companion at long wavelengths. However, many of the sec- 
ondaries are either faint, or unseen, in the SDSS spectra. In 
FigureQjJ we present the spectra of two systems where there 
is evidence for a faint companion in the spectra, and in Figure 
IT4~1 two spectra where the companion is not seen. One method 
used to determine the type of companion is to fit the spectrum 
with a two-component model comprised of a white dwarf 
plus an M-dwarf companion (Silvestri et al. 2006; Rebassa- 
Mansergas et al. 2007). This method works well when both 
components have high S/N. However, systems with very faint 
companions cannot be well decoupled when they are poorly 
constrained by the spectra. In particular, the spectra of M- 
dwarfs exhibit a degeneracy between metallicity and effective 
temperature (Bean, Fritz Benedict, & Endl 2006). This makes 
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it difficult to solve for properties of individual components of 
such systems uniquely. 

An alternative approach to constraining the secondaries is 
to use near-IR or IR data. At long wavelengths cool stars out- 
shine white dwarfs by a large factor. Therefore, to determine 
the spectral types of the secondary companions we use the 
distances derived from white dwarf models along with J,H,K 
near-IR photometry (from UKIDSS and 2MASS), and the ab- 
solute magnitudes for dwarf stars from Hoard et al. (2007). 
We matched the white dwarf locations with 2MASS detec- 
tions to well below the recommended la all-sky catalog lim- 
its (J,H,K S = 15.8, 15.1, 14.3). The bulk of our 2MASS 
matches are fainter than the formal completeness limits and 
thus should be regarded with caution. Indeed, Farihi (2009) 
noted that many white dwarf stars selected to have a near-IR 
excess from 2MASS data were found to be spurious in more 
accurate IRTF data. We also matched the white dwarfs with 
UKIDSS Large Area Survey, which has limits that are ap- 
proximately four magnitudes fainter than 2MASS. These data 
are given in Table 3. In cases where there is a near-IR detec- 
tion, the spectral type given is the average determined from 
the J,H,K measurements. If there is no detection, because 
the secondary plus white dwarfs are fainter than the limit of 
2MASS or UKIDSS, the latest type dwarf above the 2MASS 
detection threshold is given as an upper limit. 

3.4. Eclipses 

Eclipsing binary systems are very useful for determining 
properties of the associated stars. For instance, it is possi- 
ble to determine the size of the secondary using the primary 
mass, size, orbital period and the eclipse time. The depth and 
shape of an eclipse can be used to constrain the luminosity 
of the eclipsing object as well as the orbital inclination, ec- 
centricity and orientation of a binary system. White dwarf 
eclipsing systems where an object is measured photometri- 
cally during mid-eclipse will mainly be due to stellar secon- 
daries. Most eclipsing white dwarf systems discovered here 
are not detected during eclipse. 

Recently, Drake et al (2009) reported on the discovery 
of eclipsing white dwarf systems in the Catalina Real-time 
Transient Survey. The eclipsing stars were evidently late-M 
dwarfs. However, their eclipses were seen in two of four con- 
secutive images separated by 10 minutes. This fact implies 
that the eclipses lasted between 10 and 30 minutes. Based on 
Figure|2] this eclipse length suggests that the objects are larger 
than planets. Secondly, the depth of the eclipses was approx- 
imately a couple of magnitudes, suggesting that the eclips- 
ing object was fainter than the white dwarfs, but not greatly 
so. Lastly, the detection of these objects in 2MASS data con- 
firmed the companions were M dwarfs. Thus, even without 
spectra one can deduce these are WD-MD systems. 

To take advantage of the eclipsing information for each ob- 
ject, we used the light curves to determine the full length of 
the eclipse from first contact to fourth contact. We then in- 
terpolated an updated version of the Bergeron, Wesemael, & 
Beauchamp (1995) white dwarf grids to determine the mass 
and radius from the T e g and log(g) derived from the spectra. 
With this information we constrain the mass and radius of the 
secondary. That is to say, for a central transiting eclipsing 
system we have: 



where R sec , /? w( j, M wc j and M sec are the radii and masses of the 
components, Af is the eclipse time and P is the orbital period. 

Using Equation 1, we can plot the secondary star's radius 
as a function of its mass for the given white dwarf mass. By 
assuming a mass-radius relationship for companions ranging 
from planetary companions to M-dwarfs, we found the sec- 
ondary's radius by determining the intercept with equation 
1. For objects more massive than 0.1 M©, we assumed the 
mass-radius relation given by Caillautt & Patterson (1990). 
For objects less massive than 0.01 M Q we used the relation of 
Hartman et al. (2009). We note that Hartman et al.'s relation 
appears to agree well with results from planetary transit sur- 
veys where both radii and masses are known. However, in the 
brown dwarf region, between late-M dwarfs and hot jupiters, 
radii are not well understood (Burrows et al. 2001). For ob- 
jects in this region, we assume a radius of 1.5Rj, consistent 
with both late-M dwarfs and hot jupiters. In summary, we 
have: 
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(1) 



(2) 

To take into account the uncertainty in the transit times de- 
rived from the CSS light curves, we use the the longest and 
shortest possible eclipse times to determine the uncertainty in 
radius for each, combining this with the above relations. In 
addition to the uncertainty in eclipse length, there is uncer- 
tainty in the inclination of the eclipsing systems. We incorpo- 
rate this effect by determining the size of an eclipsing com- 
panion that would cause a total eclipse of the length observed 
near the companion's limb and for a central eclipse. 

In cases where a long eclipse is observed, the eclipsing ob- 
ject could be very large and simply eclipsing the white dwarf 
near its limb. However, the radius is further constrained based 
on the mass of the secondary and the brightness of the ob- 
ject in the near-IR data. For short eclipse times, the size of 
the secondary is strongly constrained by the time required for 
the secondary to completely cross the white dwarf. The radii 
derived from these data are presented in Table 3 (for systems 
with SDSS spectra and known periods). However, we caution 
that the derived companion size depends on the radius derived 
from models for the white dwarf. Comparison of models with 
values derived for Hipparcos by Provencal & Shipman (1998) 
and those presented in Pyrzas et al. (2009) suggests that the 
Bergeron et al. (1995) DA white dwarf model radii are sig- 
nificantly larger than those observed for a range of mass and 
temperature. We note that there can be considerable variation 
in white dwarf model radii depending on internal composi- 
tion (Panei et al. 2000), which in turn depends on the rate 
of the (uncertain) C 12 (a,7)0 16 reaction rate (e.g., Metcalfe, 
Salaris, & Winget 2002; see also Hansen & Liebert 2003). For 
long eclipse times, since the companion is much larger, the 
white dwarf's actual size has little effect. However, for short 
eclipses, a much smaller white dwarf leads to a significantly 
larger derived companion size. Future detailed photometric 
studies of the eclipse shape, along with radial velocities, will 
be used to constrain the inclinations and eclipse times, and 
thus masses and radii of the components in the systems, much 
more accurately than given in Table 3. 

3.5. Notes on Individual Objects 
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WD 1333+005 (LP 618-14) was first discovered as part 
of the new Luyten two-tenths catalog, as noted by Farihi et 
al. (2005a). Based on a low S/N spectrum, Farihi et al. identi- 
fied TiO bands and determined this to have a d4M companion. 
However, Farihi et al. found the white dwarf to be at 86.3 pc, 
but could not place firm limits on its temperature or distance. 
Kilic et al. (2006) independently found the white dwarf and 
determined it to have a temperature of 6100 K and a distance 
of 41 pc. They also noted the presence of the M-dwarf com- 
panion. Farihi, Hoard, & Wachter (2006) undertook HST pho- 
tometry, but did not find a resolved companion for this object. 
From their observations, they put an upper limit on the pro- 
jected orbital separation of 2.2 AU and found the white dwarf 
to be at 86 pc while the companion was found to best match 
a dM4.5 at a distance of 147 pc. The CSS data show the first 
observed eclipses in this system, and the short period (0.12 
days) confirms the companion is within 0. 1 AU of the white 
dwarf. 

SDSSJ0759, SDSSJ1502, SW Sex (PG19 12-029) and 
GK Vir (WD 1413+015) are all known cataclysmic variables 
with prominent emission in their spectra (Szkody et al. 2006). 
Of these, only SDSSJ0759 was not previously known to be an 
eclipsing system. Additionally, although CSS06833 has no 
spectrum, the variation in the CSS light curve, SDSS photom- 
etry (taken on two epochs) and 2MASS near-IR magnitudes 
suggest that this is also a CV. 

CSS080408 and CSS080502 were discovered during their 
eclipses by the Catalina Real-time Transient Survey (Drake 
et al. 2009). T e s and log(g) values of the objects were also 
derived by Pyrzas et al. (2009). For CSS080408, these val- 
ues agree within la with ours. For CSS080502, the values 
of J e ff are within la. However, the log(g) values vary by 3o\ 
CSS080502 is given to have an M3 dwarf companion by Sil- 
vestri et al. (2006), in good agreement with the M4 type de- 
rived in Table 3. 

SDSSJ0110 and SDSSJ1548 were first discovered to have 
M-dwarf companions by Silvestri et al. (2006). They derived 
M2 and M4 types for these, respectively. These two systems 
were found to be eclipsing by Pyrzas et al. (2009), who de- 
rived M4±l and M6±0.5 types. In our independent discovery 
and less detailed analysis we found M5 and >M7 types, re- 
spectively - in relatively fair agreement with Pyrzas et al.'s 
analysis. SDSSJ0110 was spectroscopically observed twice 
by SDSS. The derived T e ff and log(g) values derived from fits 
to both vary by more than the formal fit errors, but differ by 
less than 1000 K and 0.21 dex in log(g). These values are 
close to those derived by Pyrzas et al. (2009). 

SDSSJ1212 was discovered to have an M-dwarf compan- 
ion of type M2 by Silvestri et al. (2006). This system was 
found to be eclipsing by Nebot Gomez-Moran et al. (2009), 
who derived a companion of type M4±l. In our indepen- 
dent discovery and analysis we derived a companion type of 
M5. The temperature derived by autofit is consistent within 
lcr with that given by Nebot Gomez-Moran et al. (2009). 
However, the autofit log(g) value is somewhat lower (7.16 vs. 
7.6). 

The spectra of CSS09797, CSS40190, CSS21357, 
CSS38094 and CSS07125 all show clear signs of a dwarf star 
companion in their spectra. CSS07125 was found to have a 
companion of type M3 derived from the spectrum by Silvestri 
et al. (2006). This is in fair agreement with the M4 type given 
in Table 3. The companion to CSS21357 is particularly bright, 
as expected from its 2MASS near-IR magnitudes - which are 
brighter than the SDSS optical magnitudes. CSS40190 was 



photometrically observed twice by SDSS. One of these scans 
shows the object was observed in eclipse during the first three 
images in the sequence r,i,u (5r = 3,Si= 1.3 and Su = 0.3) and 
outside eclipse when the z and g images were taken. 

CSS25601 has accurate near-IR photometry from the 
UKIDSS survey that constrains the secondary to M5 type and 
exhibits emission at wavelengths > 700 nm in the SDSS spec- 
trum. However, the radius derived for this object is much 
smaller than expected for an M5 dwarf, at only 0.068 Rq. 
Although the eclipse appears to be at least three magnitudes 
deep, it is possible that this is not a complete eclipse. There- 
fore, the true size of the companion could be slightly larger 
than estimated. 

As mentioned earlier, CSS09704 lies within the SDSS 
stripe82 area that was observed many times during the SDSS- 
II supernova survey. It therefore has many multi-colour ob- 
servations. The combined CSS and SDSS light curve was 
presented in Figure [12] As the eclipse is more than two mag- 
nitudes deep in all SDSS images the companion is signifi- 
cantly fainter than the white dwarf even at red wavelengths. 
Although the eclipse time is short, the presence of a signifi- 
cant reflection effect suggests a moderately large companion. 
Therefore, this may be an eclipse near the limb of a moder- 
ately large M-dwarf. 

CSS03170, CSS41177 and CSS06653 do not show clear 
signs of M-dwarf companions in their spectra and all have 
small radii calculated by eclipse parameters and fits to SDSS 
spectra. CSS06653's companion is strongly constrained to an 
M6 type by UKIDSS near-IR photometry. The secondary to 
CSS41 177 is weakly constrained to an M6 type by a marginal 
2MASS detection and by a very small companion radius. The 
eclipse of this white dwarf is only 0.5 magnitudes, which is 
consistent with a partial eclipse of the white dwarf. However, 
the light curve is not well sampled. The transit of a small 
faint object could also produce such a signal. The compan- 
ion to CSS03170 is constrained to types later than M3 by the 
non-detection in 2MASS data. CSS03170 was observed spec- 
troscopically twice by the SDSS as shown in Figure [15] The 
temperatures and log(g) values derived from these spectra are 
consistent within 3a. 

4. DISCUSSION 

White dwarf stars are known to exist in a number of types 
of binary systems. Post-common-envelope binaries (PCEBs) 
consist of a white dwarf and main-sequence star that have 
evolved through giant branch stages of evolution and remain 
bound. Such systems are thought to evolve into CVs, where 
the matter is accreted from the main-sequence star into the 
white dwarf or an accretion disk surrounding the white dwarf. 
Such systems are relatively common and exhibit significant 
variability through outbursts (Szkody et al. 2002). A number 
of eclipsing CV systems are known to exist (see, e.g., Szkody 
et al. 2003). In such systems large eclipses often occur when 
the accreting object eclipses the accretion hot spot. CV sys- 
tems are generally characterized by strong Balmer emission 
lines. These systems evolve to shorter orbital periods as they 
lose angular momentum. CV systems are well known to ex- 
hibit a period gap in the 2-3 h range where very few systems 
have been found. As noted earlier, only a couple of objects 
lie within the CV period gap. Evolved PCEB systems are 
believed to become low-accretion-rate CVs (Gansicke et al. 
2009). These systems may undergo a period bounce where the 
low-mass secondaries "bounce" to longer periods (Paczynski 
1981; Patterson 1998). Such systems are theorized to be very 
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common (Kolb & Baraffe 1999), but only few have recently 
been found (Littlefair et al. 2008). One possible period- 
bounce system is SDSSJ123813.73- 033933.0. This system 
has a type-L4 brown dwarf companion in a 9.3hr orbit, but 
still exhibits significant Balmer emission. If the objects dis- 
covered here are period bouncers, then their lack of emission 
is unexpected in comparison to other systems. Nevertheless, it 
is possible that some of the secondaries have evolved to such 
low mass that their accretion rate is very low, and thus they 
exhibit very little emission. 

In addition to planets, white dwarfs can have a number of 
companions, including other white dwarfs (such as seen in 
binary He-atmosphere white dwarfs), AM CVn systems, neu- 
tron stars, and black holes (Roelofs et al. 2007). Systems with 
massive binary components can be seen from the presence of 
double lines, in cases of bright companions. Additionally, the 
presence of close companions gives rise to significant radial 
velocity variations as the objects orbit. Kilic et al. (2010b) dis- 
covered evidence for four binary white dwarfs systems with 
very large (> 500km/s) radial velocity variations where op- 
tical photometry was used to rule out the presence of a main 
sequence companion. Badenes et al. (2009) discovered evi- 
dence for an unseen neutron star or black hole companion to 
massive white dwarfs by following candidates selected from 
radial velocities derived from triplet spectroscopic observa- 
tions taken by SDSS. Similarly, Mullally et al. (2009) dis- 
covered two binary white dwarf systems using SDSS spectro- 
scopic data. At long orbital periods it may be possible to pro- 
duce systems where one white dwarf has cooled and causes 
significant eclipses of the other. Since most of the eclips- 
ing white dwarf systems we discovered are relatively bright, 
it is possible that some eclipsing secondaries could be much 
fainter white dwarfs. However, since the eclipses observed are 
deep, such discoveries are disfavored by the fact that the white 
dwarfs are of similar size. Even relatively low-resolution ra- 
dial velocity measurements should reveal such objects. 

Current and upcoming studies to discover and monitor ex- 
oplanets, such as Kepler and CoRoT, promise to provide 
a wealth of information about the planetary companions to 
main-sequence stars. Techniques involving radial velocities 
and planetary transits have mainly been limited to bright or 
small stars. More massive stars that have evolved to the white 
dwarf stage have not been probed for close planetary com- 
panions for three main reasons: (1) white dwarfs are mainly 
thought to have evolved through the AGB phase, and theo- 
retical studies have long suggested that close planets cannot 
survive this evolutionary phase; (2) white dwarfs are small, 
and are therefore much less likely to be seen to eclipse; (3) 
white dwarfs are generally faint and have broad absorption 
lines, which makes it difficult to obtain radial velocities with 
great precision for them. 

However, and as we have noted previously, theories regard- 
ing the survival of planetary companions during the AGB 
phase are still very uncertain, some suggesting terrestrial 
planets in our own system could survive the AGB phase, while 
others disagree. Indeed, even if planets cannot survive the 
AGB phase in close proximity to their host star, it has been 
shown that the interactions of Jupiter-like planets in planetary 
systems can give rise to unstable orbits which can later re- 
sult in close companions to white dwarfs. In addition, sdB 
stars are believed to bypass the AGB stage altogether, sim- 
ilarly to He white dwarfs. Furthermore, there is significant 
and growing evidence for recent events in white dwarf sys- 
tems involving the destruction of low-mass objects formed 



from terrestrial-type material, as well as for planets that must 
have survived the AGB phase. 

Although white dwarfs are certainly small, it is this very 
fact that makes the detection of very small transiting or eclips- 
ing planets possible. The very large signal of eclipsing events 
compensates for the faintness of white dwarfs and makes dis- 
coveries possible even with small, inexpensive, ground-based 
telescopes. The limited precision required for such detections 
means that surveys which can quickly monitor large areas of 
the sky are ideal for such discoveries. Future large transient 
surveys, such as LSST and PanSTARRS, will cover the en- 
tire sky dozens of times per year in multiple filters reaching 
stars to r = 24.5 (Ivezic et al. 2008). Such surveys are ideal 
for the discovery of large numbers of faint white dwarf can- 
didates via their distinctive colours. The synoptic nature of 
these surveys can reveal the presence of very faint and dis- 
tant white dwarfs being eclipsed or transited by exoplanets, 
or even exo-asteroids. 

The white dwarf systems that have already been observed 
in eclipsing systems, such as those presented here, may be an 
ideal place to look for transiting planets. If additional com- 
panions exist, they are likely to be close to coplanar. 

The eclipsing white dwarf systems discovered here are all 
moderately distant (200-1000 kpc), with the exception of 
WD1333+005. Many of the eclipsing white dwarf systems 
discovered are at high Galactic latitudes because most of 
the SDSS DR6 photometry concentrates on latitudes > 20°. 
These objects have projected distances > 400 pc from the 
Galactic plane, suggesting that they are part of the Galac- 
tic thick-disk population. Being in the thick disk suggests 
they are old stars, and also have average metallicity below 
solar. Additionally, the white dwarfs discovered here have 
low masses, and this suggests that their precursors had low 
masses (Weidemann 2000). Thus, this is consistent with the 
progenitors being long-lived stars. Indeed, most of the eclips- 
ing white dwarfs have relatively high temperatures. This 
may be because they have had less time to cool than more 
massive white dwarfs. The masses of these eclipsing white 
dwarfs are in fact among the lowest found among those re- 
leased in the SDSS-DR1 sample by Hu et al. (2007). How- 
ever, the masses are consistent with many determined by 
Ciechanowska et al. (2007) based on the much larger SDSS- 
DR4 white dwarf sample. Origin within the thick disk is not 
completely unexpected as Napiwotzki (2009) found that one 
third of the local population of white dwarfs comes from the 
Galactic thick disk. 

Based on CSS lightcurves, a number of eclipsing white 
dwarf systems that we have found have small companions. In 
many cases the companions are clearly M-dwarf secondaries 
because of spectroscopic features in the SDSS data and their 
measured near-IR brightness. A small number of systems ex- 
hibit no secondary in either spectrum or near-IR data. These 
companions could well be brown dwarfs or planets. However, 
further observations of these systems are required to better 
constrain the eclipse depth, shape and length. Such data will 
constrain the luminosity of the secondaries and the orbital pa- 
rameters of the systems. Additionally, deep near-IR photom- 
etry, such as available from future UKIDSS releases, can be 
used to constrain the presence of brown dwarf companions, 
and Spitzer IR observations can reach still fainter compan- 
ions. 

The fraction of white dwarfs with stellar companions is not 
well constrained by current large surveys such as the SDSS 
because they have often been serendipitously discovered in 
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searches designed to discover QSOs (Kleinman et al. 2004). 
Nevertheless, large numbers of such systems have been dis- 
covered (Silverstri et al. 2006). Under the assumption that the 
number of white dwarfs with Jupiter-size planetary compan- 
ions is the same as the number with M-dwarf companions, due 
to their relative sizes, one expects to observe eclipses for 2 to 5 
times more white dwarfs with M-dwarf companions than with 
Jupiter-size ones. Detailed follow-up of the current systems 
can constrain the number of Jupiter-size companions. How- 
ever, under the same assumption for Earth-size companions, 
we expect to find 25 to 50 times more systems with M-dwarfs, 
due to the size difference. Clearly, a few hundred systems are 
required to put significant limits on such white dwarf compan- 
ions. Given the small probability of observing eclipses, one 
must monitor many tens of thousands of white dwarfs. LSST 
should provide sufficient numbers of white dwarfs for such 
investigations. 

The CSS survey is funded by the National Aeronautics and 
Space Administration under Grant No. NNG05GF22G issued 
through the Science Mission Directorate Near-Earth Objects 
Observations Program. The CRTS survey is supported by the 
U.S. National Science Foundation under grants AST-0407448 
and AST-0909182. A.D. was supported by the National Sci- 
ence Foundation under Grant No. CNS-0540369. Support for 
M.C. is provided by MIDEPLAN's Programa Inicativa Cien- 
tifica Milenio through grant P07-021-F, awarded to The Milky 
Way Millennium Nucleus, by Proyecto Basal PFB-06/2007, 
by FONDAP Centro de Astroffsica 15010003, and by a John 
Simon Guggenheim Memorial Foundation Fellowship. This 
research has made use of the SIMBAD database, operated 
at CDS, Strasbourg, France. This research has made use of 
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FIG . 1 . — Probability of an eclipse being observable along our line-of-sight for a white dwarf with a radius of 0. 1 eclipsed by objects in the planetary range 
as a function of orbital separation. 
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FIG. 2. — Periods (P) and eclipsing times (Te) as function of orbital radius a for companions of various sizes around a white dwarf primary. The solid line 
presents the orbital period for a planet with mass M ffi orbiting a white dwarf with mass 0.6M Q . The dashed lines present the eclipse times for systems with the 
following parameters: 0.01 Rq, lMg; 0.05 Rq, 125 M©; 0.1 R0, 1000M®; 0.25 R©, O.2M0. The vertical line shows the approximate Roche stability limit. 




FIG. 3. — Selection of white dwarf candidates. Dots, SDSS DR6 sources falling within the magnitude limits noted in the text; crosses, spectroscopically 
confirmed white dwarfs from SDSS DR4 (Eisenstein et al. 2006), limited to 3000 for clarity. 
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FIG. 4. — The equatorial coordinates of the white dwarf candidates analysed in this survey, dominated by objects selected by colour from SDSS DR6 data. 
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FIG. 5. — Colours of the eclipsing white dwarfs discovered. Squares, eclipsing white dwarfs; dots, DR6 colour selected candidate white dwarfs. The dashed 
line shows the limits of the colour selection. 
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FIG. 7. — Colour distributions of white dwarfs, WD-MD binaries and CVs. The black points are SDSS DR4 white dwarfs from Eisenstein et al. (2006). Red 
squares are eclipsing white dwarfs. Blue circles are Szkody et al. (2002, 2003, 2004, 2005, 2006, 2007) cataclysmic variables. Green crosses are WD-MD 
binaries from Silvestri et al. (2006). 
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FIG. 8. — Colours of white dwarfs, WD-MD binaries and CVs. The black points are SDSS DR4 white dwarfs from Eisenstein et al. (2006). Red squares are 
eclipsing white dwarfs. Blue circles are Szkody et al. (2002-2007) cataclysmic variables. Green crosses are WD-MD binaries from Silvestri et al. (2006). 




Fig. 9.— 



QSOs in the colour-selected WD sample. The red squares are eclipsing white dwarfs. Blue triangles are variable objects in the selection region. 
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T E (mins) 

FIG. 10. — Orbital periods and eclipse times for eclipsing white dwarfs discovered herein. The solid lines shows the eclipse time (Te) for a 1.1 Mq white dwarf, 
while the dashed lines present times for 0.3 Mq white dwarfs. Typical white dwarfs will have eclipse times between the solid and dashed lines. The dotted lines 
show the bounds of the 2-3 h CV period gap. 
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FIG. 1 1. — The folded light curves of two eclipsing white dwarfs. Dots show CSS measurements, while triangles show detection upper limits. Left: white 
dwarf CSS41631. Right: white dwarf CSS06653. 
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FIG. 12. — The folded light curves of an eclipsing white dwarf, CSS09704, covered multiple times by SDSS-II supernova survey of stripe 82. Filled circles 
denote CSS unfiltered measurements, while large filled triangles show CSS upper limits. Crosses, small triangles, open boxes, open triangles, and filled boxes 
show SDSS u,g,r,i, and z measurements, respectively. 




FIG. 13. — The SDSS spectra of two eclipsing white dwarfs discovered in this study. Both spectra appear to show some evidence for a luminous companion at 
wavelengths beyond 7000 A. Left: the spectrum of CSS25601. Right: the spectrum of CSS06653. 
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FIG. 15. — The two SDSS spectra of CSS03170, a white dwarf discovered in this study. Both spectra show little evidence for a luminous companion beyond 
7000 A. 
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NOTE. — Col. (1), Object ID. Cols. (2) & (3), Right Ascension and Declination (J2000). Cols. (4)-(8), SDSS optical Gunn magnitudes. Cols. (9)-(l 1), 2MASS near-IR magnitudes 
(except where noted). Col. (12), USNO-B1.0 proper motion (Monet el al. 2003). Uncertainties in SDSS measurements are of order 0.01 magnitudes. 

t Near-IR magnitudes from UKIDSS rather than 2MASS. 
* Known CV system. 
b Drake et al. (2009). 

c Independently discovered by Pyrzas et al. (2009). 
^ Independently discovered by Nebot Gomez-Moran et al. (2009). 
e Spectroscopically confirmed with the SMARTs 1.5m telescope. 
f Not spectroscopically confirmed. 
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TABLE 2 
Eclipse Details 



ID 


HMJD E 


Period 


E L 




W L) Ljjj+K)\jj 


J3404. jyuj^o; 


A 1110?C7/1^\ 

0. Iziyjo /4(J ) 


U.U4/J 


U.UUzj 


rN*\41 177 

V.iJO'T 11// 


JjH/U.ZUjj^I / ) 


1 1 601 S4.C1 ^ 


015 


0.005 


cns91616 


S3470 ^tfYWIQ^ 

J Jt / U.jUjlH 1 y ) 


1 949S8XCM 


04 


0.002 


o VV ij CA 


S^4QC* 91 91 CI ^ 


1 3493830C9S'! 




07 


LjOi- JUU 1 


^3466 361 1 CI \ \ 


9978569C9*! 


0.0225 


0025 


rN^09797 


671 9H4^ 


9991 909C9^ 
\j.£.z.y i z,uz^z ^ 


0.0475 


0175 


SDSST0759 


JJtU7. 1 UjU^U y 


1 309337C5^ 

\). 1 JV/7JJ / \ OJ 


0.1 


0.02 


CSS401 90 


53469 9197(13^ 


1 301 1 995C40^ 

yj. 1JU1 1 Z.Z.^7(^T\7 ) 


0.045 


0.015 


CNS41631 


S3470 9634C1 5^ 


1 50X7065C1 S'i 


0.0575 


0.0025 


CSS38094 


53495.4534(33) 


0.3309896(2) 


0.033 


0.003 


CSS080408 


53470.3979(26) 


0.3820040(15) 


0.030 


0.005 


SDSSJ1502 


53479.3817(8) 


0.05890943(9) 


0.1 


0.01 


CSS06833 


53480.4230(8) 


0.16177052(8) 


0.09 


0.005 1 


CSS03170 


53464.2198(9) 


0.06509654(3) 


0.05 


0.02 


CSS09704 


53507.4564(11) 


0.1565057(2) 


0.035 


0.005 


CSS06653 


53466.3168(4) 


0.08096622(2) 


0.067 


0.002 


SDSSJ1212 


53494.6001(50) 


0.335871(3) 


0.0475 


0.0125 


CSS07125 


53464.4873(36) 


0.363497(25) 


0.03 


0.0075 


GKVir 


53464.4787(34) 


0.3443305(11) 


0.0265 


0.0085 


SDSSJ1548 


53526.2843(27) 


0.1855162(15) 


0.0375 


0.0225 


CSS080502 


53466.3302(14) 


0.1494385(25) 


0.09 


0.005 



NOTE. — Col. (1), Object ID. Col. (2), Mid-eclipse time. Col. (3), Eclipse period in days. Col. (4), Eclipse phase length. Col. (5), Uncertainty in eclipse phase length, 
t Not a confirmed white dwarf. 
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TABLE 3 
Binary System Parameters 



ID 


T e ff 


log(g) 


M w d 


D w d 




Secondary 


R-sec 




(K) 




(M ) 


(PC) 


(R©) 




(R©) 


CSS25601 


21168(240) 


7.39(0.041) 


0.37 


452 


0.020 


M5 


0.05 - 0.068 


LjjU? ly I 


1 Q 1 AQC\ (*(W 
ly my\ 1DU ) 


7 iivn m 1 ■! 


U. J4 


ion 


n mo 


> M J 


U.UoZ - U.J 






7 7/1/0 (MQ\ 


U.JU 


1 06 


n n i fi 

U.U1U 


M7 


U.U 1 - U.J 


CSS21357 


16092(195) 


7 27C0 035"! 


0.32 


287 


0.022 


M4 




CSS21055 


14462(170) 


7.66(0.038) 


0.49 


242 


0.017 


> M6 




CSS41631 


24526(274) 


7.27(0.040) 


0.35 


497 


0.023 


M5 


0.124-0.3 


CSS41177 


20930(156) 


7.37(0.027) 


0.36 


326 


0.021 


M6 


< 0.025 


CSS080408 


32041(224) 


7.20(0.051) 


0.36 


663 


0.025 


M3 


0.1-0.35* 


CSS38094 


28442(209) 


7.58(0.036) 


0.49 


467 


0.019 


M5 


0.018-0.088* 


CSS03170 


34951(412) 


7.10(0.079) 


0.34 


925 


0.027 


>M3 


CSS03170 


39040(302) 


7.30(0.029) 


0.40 


876 


0.024 


>M3 


0.018-0.088* 


CSS080502 


17120(163) 


6.96(0.036) 


0.26 


316 


0.026 


M4 


0.21-0.5* 


SDSSJ0110 


24621(180) 


7.40(0.023) 


0.38 


262 


0.020 


M5 


0.41 -0.5t.t 


SDSSJ0110 


23841(180) 


7.29(0.026) 


0.35 


274 


0.022 


M5 


0.41 -0.5*' 1 


CSS06653 


13589(114) 


7.39(0.025) 


0.34 


213 


0.020 


M6 


0.094-0.128 


SDSSJ1212 


17543(92) 


7.16(0.022) 


0.30 


264 


0.024 


M6 


0.126-0.5* 


CSS07125 


29268(510) 


7.34(0.093) 


0.38 


823 


0.022 


M4 


0.088 - 0.5 


SDSSJ1548 


9698(69) 


7.77(0.088) 


0.43 


193 


0.018 


>M7 


0.032 - 0.5* 


CSS09704 


29302(392) 


7.30(0.081) 


0.37 


744 


0.022 


>M4 


0.06 - 0.085 



NOTE. — Col. (1), Object ID. Col. (2), white dwarf effective temperature and formal fit error. Col. (3), surface gravity and fit error. Col. (4), white dwarf mass. Col. (5), distance 
derived from WD models. Col. (6), white dwarf radius based on models. Col. (7), secondary MK type based on distance from near-IR data and models. Col. (8), secondary radius. 

T Indepentently discovered by Pyrzas et al. (2009) and Nebot Gomez-Moran et al. (2009). 

* From Drake et al. (2009). 

* Object has two spectra. Results were derived separately for each. 



